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Abstract 

Recently, several groups identified a tentative 7-ray line signal with energy 
~ 130 GeV in the central Galaxy from the Fermi-LAT data. Such a 7-ray line 
can be interpreted as the signal of dark matter annihilation. However, the 
offset ~ 220 pc (1.5°) of the center of the most prominent signal region from 
the Galactic center Sgr A* has been thought to challenge the dark matter 
annihilation interpretation. Considering the fact that such a 130 GeV 7-ray 
line signal consists of only ~ 14 photons, we suggest that the "imperfect" 
consistency of these photons with the expected dark matter distribution is 
due to the limited statistics. The offset will be smaller as more signal photons 
have been collected in the near future. Our Monte Carlo simulation supports 
the above speculation. 

High energy 7-ray line is of extreme interest in search for the signal of 
dark matter (DM) annihilation or decay. Recently, via analyzing the publicly 
available Fermi-LAT 7-ray data, Bringmann et al. [H] and Weniger j2| found 
weak evidence for a monochromatic 7-ray line with energy ~ 130 GeV. Later 
independent analyses carried out by a few groups confirmed the existence 
of the 130 GeV 7-ray excess, and the signal has been found to be even 
more prominent 0, 0, [sj. This result can be interpreted by ~ 130 GeV DM 
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annihilation, with annihilation cross section {av)^^^^^ ~ 10 
a cuspy density profile such as Navarro-Frenk- White (NFW, 
0. 

Much attention was paid on this line signal in the community. Many 
models were proposed to explain this line structure, either by DM [sj or 
astrophysical sources jo], [loj. It was also suggested to constrain the DM 



scenarios with the continuum 7-rays or antiprotons or to test the line 



postulation with high energy resolution detectors 12j . Moreover, the spectra 



of the sum of cosmic ray electrons and positrons detected by ATIC and 
PAMELA both showed small wiggle-like structures at ~ 100 GeV 0, Q, 
which could be the result of the annihilation of ~ 140 GeV DM particles into 



electrons/positrons, in accordance with the 130 GeV 7-ray line 15 



Among current relevant data analysis works, the morphology of the po- 
tential line signal is still in debate 0, 0, Isl- It is very attractive that Su 
& Finkbeiner identified that the signal region lies basically in the Galactic 
center region, and the detection significance is quite high (exceeds 5a) jsj. 
The former character is just expected in the DM scenario, while a confidence 
level > 5a, if confirmed, is encouraging to approach a discovery/detection. 
The problem is however that the signal region has a center deviating from 
the Galactic centre (Sgr A*) considerably by a distance ~ 220 pc (or angle 
~ 1.5°)|^. The result seems to be at odds with the DM models in which the 
signal region is expected to be centered at (£, b) = (0°, 0°), where (£, b) are 
the Galactic longitude and latitude, respectively. Indeed, such a puzzle has 
been thought to be one of the strongest arguments against the DM origin 



of the 7-ray line signal |10l . Il6| . In this Letter, considering the fact that the 
current 130 GeV 7-ray line signal consists of only ~ 14 photons, we try to 
provide a statistical interpretation of the spatial distribution. 

For such a purpose we carry out the Monte Carlo simulation of the ar- 
rival direction of the photons produced by the annihilation of DM particles. 
Following jsf we adopt the Einasto DM density profile in this work 



p(r) = ps exp 



a 



(1) 



where a = 0.17, 



20 kpc and ps ~ 0.06 GeV cm 



^For the signal region with highest significance (~ 4.5(t, i.e., their central region) iden- 
tified in [s^l, an offset ~ 1.2° was reported, too. 
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Figure 1: The inclination-angle averaged point spread function of Fermi-LAT at 130 GeV. 

The possibility of detecting one photon at the location {£, b) is propor- 
tional to the J-factor 

J<x J dsp'iris)), (2) 

where r = (s^ + — 2srQ cos £ cos b)^^"^ is the Galactocentric distance, Tq ~ 
8.5 kpc is the distance from the Sun to the Galactic center and s is the line 
of sight distance. 

In reality the "observed" (reconstructed) direction of the photons will 
deviate from the "real" location, and the deviation is a function of an incident 
photon's energy and inclination angle. Such an effect is known as the point- 
spread function (PSF) of the instrument. For Fermi-LAT, the PSF function 
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Figure 2: Left: scattering plot of (.^o, ^o) in the 10^ simulations; right: scattering plot of 
offset angle r versus the elongation rate aijuh- 

can be found on the websit^. In this work, we fix the photon energy to be 
130 GeV and average the PSF in different inclination angles. The derived 
PSF function j[ff) is shown in Fig. [H where the normalization 27r J f{0)d6 = 
1 has been adopted. For each photon in the simulation, we re-generate the 
"observed" direction which may deviate 6 from the "real" direction with 
probability f{9). 

We simulate 10000 observations with = 14 photons each. These pho- 
tons are assumed to come from an angle ^ < 5° around the Galactic center 
since all the 14 signal photons identified in~[5[ were from such a very com- 
pact region, where cos^ = cos £ cos 6 (for the Einasto DM density profile, just 
about 1/4 of the signal photons will be from ^ < 5°, implying that there are 
many more signal photons from larger angles, which however may have been 
hidden behind the dense background). The average center of the photons is 
estimated to be 

N N 

io = Y,ii/N, bo = Y,k/N. (3) 

i=l 1=1 

The distribution of the resulting (£o, &o) is presented in the left panel of Fig. 

m 

The offset of the morphology center compared with the Galactic center 
is simply r = a/^q + ^o- We also investigate the asymmetric property of the 



■^http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/docunientation/Cicerone/Cicerone_LAT_IRFs/IRF_PSF. html 



photon map, through defining the elongation rate 
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In the right panel of Fig. |2]we show the distribution of r and otjoh- Lines in 
this figure present the la and la contours. It is shown that an offset of about 
1.5° revealed by the data is consistent with the canonical DM distribution 
within 2(7 confidence level. Specifically the probability of r > 1.5° is about 
2% (or 2.3cr), for the case iV = 14. Our prediction will be directly tested by 
the ongoing and upcoming high energy observations. The results also show 
that a potential asymmetry between the longitude and latitude directions Q 
may also appear due to the limited statistics. 

With the increase of photon statistics, we would expect the deviation of 
the morphology center from the real center to decrease. The resulting distri- 
butions of the parameters for N = 30, 50 and 100, still for 10000 simulations, 
are shown in Fig. [3l Given more and more photons we find that the fiuc- 
tuation of the morphology center becomes smaller and smaller, as expected. 
The probability for a large offset of the morphology center is accordingly 
much smaller. For = 30 we have P{r > 1.5°) = 0.02%. This probability 
decreases to 8 x 10~^ for = 50, and become much smaller for = 100. 
Obviously, the morphology will be more symmetric given more photons are 
detected. 

In summary we have shown that the spatial distribution of 130 GeV 7- 
ray line signal identified in [sl can be consistent with the DM annihilation 
model and the offset of the signal region from the Galactic centre is likely 
caused by the limited statistics. The upcoming high energy resolution de- 
tectors such as DArk Matter Particle Explorer (DAMPE) and CALorimetric 
Electron Telescope (CALETfl will be more powerful in identifying the line- 
like 7-ray signal [e.g., |l2|. These two detectors however have an effective 
area smaller than Fermi-LAT, and the total photons detectable should also 
be fewer. We thus do not expect to get a perfect coincidence of the signal 
region with the expected DM distribution even in future observations. In 
this Letter the distribution of dark matter particles in the central Galaxy is 
assumed to be centered at {i, h) = (0°, 0°). If it is not the case an offset 
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Figure 3: Same as Fig. H but for iV = 30 (top), TV = 50 (middle) and N = 100 (bottom) 
respectively. 
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lager than that predicted in our simulation is likely and the current spatial 
distribution of the signal photons might be better reproduced. 
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